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ABSTRACT 

We present the first comparison of virial masses of galaxy clusters with their Sunyaev-Zel'dovich 
Effect (SZE) signals. We study 15 clusters from the Hectospec Cluster Survey (HeCS) with 
MMT/Hectospec spectroscopy and published SZE signals. We measure virial masses of these clusters 
from an average of 90 member redshifts inside the radius rioo ■ The virial masses of the clusters are 
strongly correlated with their SZE signals (at the 99% confidence level using a Spearman rank-sum 
test). This correlation suggests that Ygz can be used as a measure of virial mass. Simulations predict 
a powerlaw scaling of Ysz oc Af^ 00 with a ssl.6. Observationally, we find a=l.ll±0.16, significantly 
shallower (given the formal uncertainty) than the theoretical prediction. However, the selection func- 
tion of our sample is unknown and a bias against less massive clusters cannot be excluded (such a 
selection bias could artificially flatten the slope). Moreover, our sample indicates that the relation 
between velocity dispersion (or virial mass estimate) and SZE signal has significant intrinsic scatter, 
comparable to the range of our current sample. More detailed studies of scaling relations are therefore 
needed to derive a robust determination of the relation between cluster mass and SZE. 
Subject headings: galaxies: clusters: individual — galaxies: kinematics and dynamics — cosmology: 
observations 



1. INTRODUCTION 

Clusters of galaxies are the most massive virialized 
systems in the universe. The normalization and evo- 
lution of the cluster mass function is therefore a sen- 
sitive probe o f the growth of s t ructure and thus cos- 
molog y (e.g.. iRines et all 120071 120081: iVikhlinin et al 



in nu merical simulations (lDiaferiol Il999t lEvrard et alJ 
2008), and recent results from hydrodynamical simula- 
tions indicate that virial mas ses may have scatter as 
small as -5% (|Lau et al.ll2010h . ' 
Previous studies have compared SZE signals to hydro 



2009; IHenrv et all 120091: IMantz et alJ 120081: IRozo eTal 
20081 and references therein). Many methods exist 



static X-ray masses (|Bonamente et al . 2008; Plagg e et al 
2010) and gravitational lensing masses ( jMarrone et al 



20091 hereafter M09). Here, we make the first compar- 



to estimate cluster m asses, includin g dynamical masses 
from either galaxie s ( Zwick^ 119371 ) or intracluster gas 
(e.g.. iFabricant et al.l I1980T). gravitat ional lensing (e.g., 
ISmith et al .ll2005HRic hard et alj2010l). and th e Sunvaev- 
Zel'dovich effect (SZE iSunvaev fc Zeldovichl IT972T) . In 
practice, these estimates are often made using simple 
observables, such as velocity dispersion for galaxy dy- 
namics or X-ray temperature for the intracluster gas. 
If one of these observable properties of clusters has a 
well-defined relation to the cluster mass, a large survey 
can y ield tight constraints on cosmological parameters 
(e.g., iMaiumdar fc Mohrl l200l . There is thus much 
interest in identifying cluster obs ervables that exhibit 
tight scaling rela tions with mass (jKravtsov et alJ 120061 : 
IRozo et "all 120081 ). Nu merical simulation s indicate that 
X-ra y gas observables (jNagai et alJ 120071) and SZE sig- 
nals (|Motl et al.ll2005D are both candidates for tight scal- 
ing relations. Both met hods are beginning to gain ob- 
serva t ional support (e.g. . IHenrv et alJl2009b lLopes et aD 
120091 : IMantz et all 120091 : lLocutus Huang et alJ 12009). 
Dynamical masses from galaxy velocities are unbiased 
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ison between virial masses of galaxy clusters and their 
SZE signals. We use SZE measurements from the lit- 
erature and newly-measured virial masses of 15 clus- 
ters from extensive MMT/Hectospec spectroscopy. This 
comparison tests the robustness of the SZE as a proxy 
for cluster mass and the physical relationship between 
the SZE signal and cluster mass. Large SZ cluster sur- 
veys are under way and are beginning: to yield cosmologi- 
cal constraints dCarlstrom et al.l2010rlHiiicks et alJ2010t 
iStaniszewski et al.ll2009D . 

We assume a cosmology of fi m =0.3, fiA=0.7, and 
-£/o=70 km s _1 Mpc -1 for all calculations. 

2. OBSERVATIONS 

2.1. Optical Photometry and Spectroscopy 

We are completing the Hectospec Cluster Survey 
(HeCS), a study of an X-ray flux-limited sample of 53 
galaxy clusters at moderate redshift with extensive spec- 
troscopy from MMT/Hectospec. HeCS includes all clus- 
ters with ROSAT X-ray fluxes of fx > 5 x 10" 12 erg 
s" 1 at r0.5-2.0lkeV from the Bright Clus ter Survey (BCS 
Ebeling et alTl9 98) or REFLEX survey (jBohringer et alJ 
20041 ) with optic al imaging in the Sixth Data R elease 
(DR6) of SDSS (|Adelman-McCarthv et afll2008l) . We 
use DR6 photometry to select Hectospec targets. The 
HeCS targets are all brighter than r=20.8 (SDSS cata- 
logs are 95% complete for point sources to r«22.2). Out 
of the HeCS sample, 15 clusters have published SZ mea- 
surements. 
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2.1.1. Spectroscopy: MMT/Hectospec and SDSS 

HcCS is a spectroscopic survey of clusters in the red- 
shift range 0.10<,z<0.30. We measure spect ra with 
the Hectospec instrument (Fabricant ct al. 20051) on the 
MMT 6.5m telescope. Hectospec provides simultaneous 
spectroscopy of up to 300 objects across a diameter of 
1°. This telescope and instrument combination is ideal 
for studying the virial regions and outskirts of clusters 
at these redshifts. We use the red sequence to preselect 
likely cluster members as primary targets, and we fill 
fibers with bluer targets (Rines et al. in prep, describes 
the details of target selection) . We eliminate all targets 
with existing SDSS spectroscopy from our target lists but 
include these in our final redshift catalogs. 

Of the 15 clusters studied here, one was observed with a 
single Hectospec pointing and the remaining 14 were ob- 
served with two pointings. Using multiple pointings and 
incorporating SDSS redshifts of brighter objects mitigate 
fiber collision issues. Because the galaxy targets are rel- 
atively bright (r<20.8), the spectra were obtained with 
relatively short exposure times of 3x600s to 4x900s under 
a variety of observing conditions. 

Figure Q] shows the redshifts of galaxies versus their 
projected clustrocentric radii for the 15 clusters stud- 
ied here. The infall patterns are clearl y present in al l 
clusters. We use the caustic technique (|Diaferiolll999D 
to determine cluster membership. Briefly, the caustic 
technique uses a redshift-radius diagram to isolate clus- 
ter members in phase space by using an adaptive ker- 
nel estimator to smooth out the galaxies in phase space, 
and then dete rmining the edges of this distribution (see 
IDiaferiol 120091 for a recent review). This technique has 
been successfully applied to optical studies of X-ray clus- 
ters, and yields cluster mass estimates in agreement 
with estimat es from X-ray obse r vations and gravitational 
lensing (e.g.. IRines et al]|2003t iBiviano fc Girardll 120031: 
Diafe rio et alj 120051: IRines fc Diaferiol 120061: IRines et al l 
20071 and references therein). 

We apply the prescription of IDanese et al.l ((19800 to 
determine the mean redshift cz & and projected velocity 
dispersion a p of each cluster from all galaxies within the 
caustics. We calculate a p using only the cluster members 
projected within rioo estimated from the caustic mass 
profile. 

2.2. SZE Measurements 

The SZE detect ions are primarily from 
iBonamente et al.l (|2008l hereaft er B08), supplemented 
by three measurements from iMarrone et al.l ([2009L 
hereafter M09). Most of the SZ data were obtained with 
the OVRO/BIMA arrays; the additional clusters from 
M09 were observed with the Su nyaev-Zel'dovich Array 
(SZA; e.g.. lMuchovei et alJl2051 . 

Numerical simulations indicate that the integrated 
Compton y-parameter Y sz has smaller sca tter than the 
peak y-decrement y pe ak (jMotl et al.l [ 2005f ), so B08 and 
M09 report only Ysz- Although y pea k should be nearly 
independent of redshift, Ysz depends on the angular size 
of the cluster. The quantity Y$zD\ removes this depen- 
dence. Thus, we compare our dynamical mass estimates 
to this quantity rather than y pea k or Ysz- Tabled] sum- 
marizes the SZ data and optical spectroscopy. 

It is also critical to determine the radius within which 



Ysz is determined. B08 use r25oo> the radius that en- 
closes an average density of 2500 times the critical den- 
sity at the cluster's redshift; r^soo has physical values of 
300-700 kpc for the massive clusters studied by B08 (470- 
670 kpc for the subsample studied here). M09 use a phys- 
ical radius of 350 kpc because this radius best matches 
their lensing data. 

To use both sets of data, we must estimate the con- 
version between Ys z{r250o) measured within r25oo and 
Ysz{f = 350 kpc) measured within the smaller radius 
r =350 kpc. There are 8 clusters analyzed in both B08 
and M09 (5 of which are in HeCS). We perform a least- 
squares fit to Ys z(f250o) ~ Ysz(r = 350kpc) to determine 
an approximate aperture correction for the M09 clusters. 
We list both quantities in Table 1. 

3. RESULTS 

We examine two issues: (1) the strength of the corre- 
lation between SZE signal and the dynamical mass and 
(2) the slope of the relationship between them. Figure [2] 
shows the Ysz — &p relation. Here, we compute a p for all 
galaxies inside both the caustics and the radius ri o, c de- 
fined by the caustic mass profile [r$ is the radius within 
which the enclosed density is 5 times the critical density 
Pc{z)]- 

Because we make the first comparison of dynami- 
cal properties and SZE signals, we first confirm that 
these two variables are well correlated. A nonparametric 
Spearman rank-sum test (one-tailed) rejects the hypoth- 
esis of uncorrelated data at the 98.4% confidence level. 
The strong correlation in the data suggests that both a p 
and YszDj, increase with increasing cluster mass. 

Hydrodynamic numerical simulations indicate that 
Ysz (integrated to r 50 o) scales with cluster mass as 
Ysz oc M£q , where a=1.60 with radiative cooling and 
star formation, and 1.61 for simulations with radiative 
cooling, star formation, an d AGN feedback (a— 1.70 for 
non-radiative simulations. iMotl et al.l 120051 ). Combin- 
ing this result with the virial scaling relation of dark 
matter particles, a p oc M 2 ° 336±0 003 (jEvrard et alJl2008| ). 
the expected scaling is Ysz oc er 4 - 76 (we assume that 
Mioo oc M500). The right panels of Figure [2] shows this 
predicted slope (dashed lines). 

The bisector of the least-squares fits to the data has 
a slope of 2.94 ± 0.74, significantly shallower than the 
predicted slope of 4.8. 

We recompute the velocity dispersions <j p .a for all 
galaxies within one Abell radius (2.14 Mpc) and in- 
side the caustics. Surprisingly, the correlation is slightly 
stronger (99.4% confidence level). This result supports 
the idea that velocity dispersions computed within a 
fixed physical radius retain strong correlations with other 
cluster observables, even though we measure the velocity 
dispersion inside different fractions of the virial radius 
for clusters of different masses. Beca use cluster veloc- 
ity dispersions decline with radius (e.g. IRines et al.ll2003l : 
IRines fc Diaferioll2006l ). a Pj A may be smaller than cr Pi ioo 
(measured within rioo.c) for low-mass clusters, perhaps 
exaggerating the difference in measured velocity disper- 
sions relative to the differences in virial mass (i.e., <j p .a of 
a low-mass cluster may be measured within 2rioo while 
fj Pi A of a high-mass cluster may be measured within rioo; 
the ratio a Pt A of these clusters would be exaggerated rel- 
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Fig. 1. — Redshift versus projected clustrocentric radius for the 15 HeCS clusters studied here. Clusters are ordered left-to-right and 
top-to-bottom by decreasing values of Ys z f ^(^2500) ■ The solid lines show the locations of the caustics, which we use to identify cluster 
members. The Hectospec data extend out to ~8 Mpc; the figure shows only the inner 4 Mpc to focus on the virial regions. 



ative to the ratio c^ioo)- Future cluster surveys with 
enough redshifts to estimate velocity dispersions but too 
few to perform a caustic analysis should still be sufficient 
for analyzing scaling relations. 

Because of random errors in the mass estimation, the 
virial mass and the caustic mass within a given radius 
do not necessarily coincide. Therefore, the radius rioo 
depends on the mass estimator used. Figure [5] shows 
the scaling relations for two estimated masses Mioo, c and 
Mioo,?;; Mioo.c is the mass estimated within rio 0iC (where 
both quantities are defined from the caustic mass profile) , 
and Mioo « is the ma,ss 6stinicit6cl within rioo v 

(both 

quantities are estimat ed with the virial theorem, e.g., 
IRines fc Diaferidl20 06). including galaxies projected in- 
side ?Too,u- Similar to a p , there is a clear correlation 
between M W o, v and YszD\ (99.0% confidence with a 
Spearman test). The strong correlation of dynamical 
mass with SZE also holds for Mioo, c estimated directly 



from the caustic technique (99.8% confidence). 

The bisector of the least-squares fits has a slope of 
1.11 ± 0.16, again significantly shallower than the pre- 
dicted slope of 1.6. This discrepancy has two distinct 
origins. By looking at the distribution of the SZE sig- 
nals in Figure [21 we see that, at a given velocity disper- 
sion or mass, the SZE signals have a scatter which is a 
factor of ~2. Alternatively, at fixed SZE signal, there 
is a scatter of a factor of ~2 in estimated virial mass. 
Unless the observational uncertainties are significantly 
underestimated, the data show substantial intrinsic scat- 
ter. Moreover, this scatter is comparable to the range of 
our sample and, therefore, the error on the slope derived 
from our least-sq uares fit to the data is l ikely to be un- 
derestimated (see lAndreon fc Hurnll20ldl . for a detailed 
discussion of a Bayesian approach to fitting relations with 
measurement uncertainties and intrinsic scatter in both 
quantities) . 
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TABLE 1 

HeCS Dynamical Masses and SZE Signals 
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z 


a p 
km s 
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1 
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1 

1 


A1689 


0.1844 


886 + S 


15.3±1.4 


9.44±5.66 


7.51±0.60 


1.50±0.14 


1 


A1763 


0.2315 


1042±™ 


16.9±1.6 


12.6±1.5 


3.10±0.32 


0.46±0.05 a 


2 


A1835 


0.2507 


1046+^ 


19.6±1.6 


20.6±0.3 


6.82±0.48 


1.37±0.11 


1 


A1914 


0.1659 


698l« 


6.70±0.57 


6.21±0.21 




1.08±0.09 


1 


A2111 


0.2290 




4.01±0.41 


4.77±1.23 




0.55±0.12 


1 


A2219 


0.2256 


915+g 


12.8±1.0 


12.0±4.7 


6.27±0.26 


1.19±0.05 a 


2 


A2259 


0.1606 


735±g 


5.59±0.60 


4.90±1.69 




0.27±0.10 


1 


A2261 


0.2249 


7251^ 


7.13±0.83 


5.10±2.07 
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1 


RXJ2129 


0.2338 


684±g 


4.31±0.57 


2.94±0.13 
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1 



Note. — "Extrapolated to r25oo using the best-fit relation between Yg z D^(350kpc) and Ys z D 2 A (r25Qo) f° r eight clusters in common 
between B08 and M09. 

Note. — Redshift z and velocity dispersion a v are computed for galaxies defined as members using the caustics. Masses Afioo.u and 
Mioo, c are evaluated using the virial mass profile and caustic mass profile respectively. 

Note. — REFERENCES: SZE data are from (1) Bonamente et al. 2008 and (2) Marrone et al. 2009. 



Our shallow slopes may also arise in part from the fact 
that our sample, which has been assembled from the lit- 
erature and whose selection function is difficult to deter- 
mine, is likely to be biased against clusters with small 
mass and low SZE signal. Larger samples should deter- 
mine whether unknown observational biases or issues in 
the physical understanding of the relation account for 
this discrepancy. 

4. DISCUSSION 

The strong correlation between masses from galaxy dy- 
namics and SZE signals indicates that the SZE is a rea- 
sonable proxy for cluster mass. B08 compare SZE signals 
to X-ray observables, in particular the temperature Tx of 
the intracluster medium and Y x = M aas Tx, wher e M gas 
is the mass of the ICM (see also lPlagge et all2010D . Both 
of these quantities are measured within rsoo, a signifi- 
cantly smaller radius than rioo where we measure virial 
mass. M09 compare SZE signals to masses estimated 
from gravitational lensing measurements. The lensing 
masses are measured within a radius of 350 kpc. For the 
clusters studied here, this radius is smaller than ^2500 
and much smaller than rinn. Numerical simulations indi- 
cate that the scatter in masses mea sured within an over- 
density S decreases as S decreases () White) 120021) . largely 
because variations in cluster cores are averaged out at 
larger radii. Thus, the dynamical measurement reaching 
to larger radius may provide a more robust indication 
of the relationship between the SZE measurements and 
cluster mass. 

The YszD 2 A — Mi ens data presented in M09 show a 
weaker correlation than our optical dynamical properties. 
A Spearman test rejects the hypothesis of uncorrelated 
data for the M09 data at only the 94.8% confidence level, 
compared to the 98.4-99.8% confidence levels for our op- 
tical dynamical properties. One possibility is that M\ ens 



is more strongly affected by substructure in cluster cores 
and by line-of-sight structures than are the virial masses 
and velocity dispersions we derive. 

Few measurements of SZE at large radii (> rsoo) are 
currently available. Hopefully, future SZ data will allow 
a comparison between virial mass and Ysz within similar 
apertures. 

5. CONCLUSIONS 

Our first direct comparison of virial masses, velocity 
dispersions, and SZ measurements for a sizable clus- 
ter sample demonstrates a strong correlation between 
these observables (98.4-99.8% confidence). The SZE sig- 
nal increases with cluster mass. However the slopes of 
both the Y sz - cr relation (Y sz oc a^ 9 ^ - 74 ) and the 

Ysz — Mioo relation (Ysz oc All Q Q 1±0 ' 16 ) are significantly 
shallower (given the formal uncertainties) than the slopes 
predicted by numerical simulations (4.76 and 1.60 respec- 
tively). 

This result may be partly explained by a bias against 
less massive clusters that could artificially flatten our 
measured slopes. Unfortunately, the selection function 
of our sample is unknown and we are unable to quan- 
tify the size of this effect. More importantly, our sample 
indicates that the relation between SZE and virial mass 
estimates (or velocity dispersion) has a non-negligible in- 
trinsic scatter. A complete, representative cluster sample 
is required to robustly determine the size of this scatter, 
its origin, and its possible effect on the SZE as a mass 
proxy. 

Curiously, Ysz is more strongly correlated with both 
(T p and Afioo than with Mi ens (M09). Comparison of 
lensing masses and cluster velocity dispersions (and virial 
masses) for larger, complete, objectively selected samples 
of clusters may resolve these differences. 

The full HeCS sample of 53 clusters will provide a large 
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Fig. 2. — Integrated S-Z Compton parameter Y$zD\ versus dynamical properties for 15 clusters from HeCS. Left panels: SZE data 
versus virial mass Mioo estimated from the virial mass profile (top) and the caustic mass profile (bottom). Solid and open points indicate 
SZ measurem ents from B08 an d M09 respectively. The dashed line shows the slope of the scaling predicted from numerical simulations: 
Y$z °t M 1,6 (Motl ct al. 2005), while the solid line shows the ordinary least-squares bisector. Arrows show the aperture corrections to 
the SZE measurements (see text). Right panels: SZE data versus projected velocity dispersions measured for galaxies inside the caustics 
and (top) inside rioo, c estimated from the caustic mass profile and (bottom) inside the Abell radius 2.14 Mpc. The dashed line shows the 
scaling predicted from simulations: Ysz oc M 1 - 6 ( Motl et al. 2005) and a oc M 33 IIEvrard et al.ir20"08l ). The solid line shows the ordinary 
least-squares bisector. Data points and arrows are defined as in the left panels. 



sample of clusters with robustly measured velocity dis- 
persions and virial masses as a partial foundation for 
these comparisons. 

We thank Stefano Andreon for fruitful discussions 
about fitting scaling relations with measurement errors 



and intrinsic scatter in both quantities. AD gratefully 
acknowledges partial support from INFN grant PD51. 
We thank Susan Tokarz for reducing the spectroscopic 
data and Perry Berlind and Mike Calkins for assisting 
with the observations. 
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